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Abstract 

Recently, there emerges different versions of beta function and hypergeometric functions containing extra 
parameters. Gaining enhghtenment from these ideas, we will first introduce a new extension of generalized 
hypergeometric function and then put forward some fundamental results in the paper. Next, we will de- 
rive some properties of certain functions like extended Gauss hypergeometric functions, extended Appell's 
hypergeometric functions F^^^', F2'^'', and extended Lauricella's hypergeometric function f[^'j^^j,F^'j|^^j, 
including transformation formulas, finite sum representations, Mellin-Barnes type integral representations 
and recurrence relations. Moreover, by using some new integral representation which will be presented 
in this paper, a Hardy-Hilbert type inequality involving extended Gauss hypergeometric functions will be 
established. 

Keywords: Extended beta function. Generalized hypergeometric functions, Appell's hypergeometric 
functions, Lauricella's hypergeometric functions, Ramanujan's Master Theorem, Method of brackets. 



1. Introduction and preliminaries. 

In the eighteenth century, L. Euler (1707-1783) concerned himself with the problem of interpolating 
between the numbers 



n! = / e-H^dt, n = 0,1,2,- •• 





with nonintergral values of n. This problem led Euler in 1729 to the now famous gamma function, a 
generalization of the factorial function that gives meaning to xl when x is any positive number. His result 
can be extended to certain negative numbers and even to complex numbers The integral representation 
of now widely accepted gamma function P [x] is 

r(x)=/ f-^e-^dt, $HM >0. (1) 



JO 

In 1994, by inserting a regularization factor e^*'^ , Chaudhry and Zubair Q have introduced the following 
extension of gamma function 

Tp [x] = ^ t"-^ exp (^-t - dt, D\ (b) > 0. (2) 

For (b) > this factor removes the singularity coming from the t = limit and for b = reduces to the 
original gamma function defined by (1). In 1997, Chaudhry et al. Q presented the following extension of 
Euler's beta function 



Bb(x,y) ^^'t^-i (1 -t)y-^xp dt, m[h)>0. 



(3) 



Afterwards, Ghaudhry et al. Q used Bb [x,y) to extended Gauss hypergeometric function and Kummer 
confluent hypergeometric function as follows 

Fb(a,b;c;z)= ^ (a)^ EfbT-V > 0; |z| < 1; 9^ (c) > 9^ (b) > 0) (4) 

Tri=0 ^ ' ' 
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0,(b;c;z)^X ^^Mbr-'b/^ S' (b>0;9^(c)>9^(b)>0) (5) 

m=0 ^ ' ' 

where (a)^ denotes the Pochhammer symbol defined, in terms of gamma functions, by 

_ r(a + m) _ f 1 m = 0; a e C\{0} 

r(a) ~ |a(a+1){a + 2)---(a + m-1) m e N; a e C. 

For b = 0, functions (4) and (5) reduces to the usual hypergeometric functions. 

In a similar manner, M. Ali Ozarslan and E. Ozergin [2l| defined the extension of Appell's functions 
Fi (a, b, c; d;x,y; b), F2 (a, b, c; d, e; x,y; b) and extended Lauricella's hypergeometric function f\^\ (q, b, c, d; e;x,y, z) 

In 2011, E. Ozergin, M. Ali Ozarslan, and A. Altin [22,] introducing the following generalizations 

r^*'P' [x) =^ t''-^Fl (^a;|3;-t-^^ dt [m[cc) > 0, $K((3) >0,in(p) >0,5H(x] > 0] (6) 

B^^' (x,y) t"-! (1 -t)^-^ iF, (^a;(3;^^^) dt. (7) 

[m[a)>0, fR (P) > 0, fR (p) > 0, 9^ (x) > 0, 5n (y ) > 0) 

In this paper, our extension mainly based on using the following generalization of gamma and beta 
functions. 

Definition 1.1. fl^ . p. 24-3]: Let a function & ({Kij^gp,^ ;z) be analytic within the disk |z| < R (0 < R < 00) 
and let its Taylor- Maclaurin coeficients be explicitly deonted by the sequence . Suppose also that the 

function @ ({Kij^^pj^ ;z) can be continued analytically in the right half-plane ^ft (z] > with the asymptotic 
property given as follows: 

e(K^;z)^e({Kj,,, ;z) = |^^-°^^^ , (|z|<R;0<R<oo;ko = 1) 

VI LJieNo' ; iMoz'^exp^z) [1 +0(1)] (9^ (z) ^ 00; Mq > 0; o) e C) ^' 

for some suitable constants Mq and cu depending essentially on the sequence [Kl\^^^^ . We can define extended 
Gamma function ^ (z) and the extended Beta function B^*^^^ (a, (3;b) by 

tI^'\z]^ j\--'@[[Ky}--X-'^di\ (z) > 0; $n (b) > 0) (9) 



and 



4^^' (a, (3;b) = j'^ t«-i (1 -tf-' e (^M;-;^^^) dt 



(10) 



(min{<H(a),<H((3)}>0; fR(b)>0) 
By introducing one additional parameter d with $H(d) > 0, we have 

-1 



^dt'^-P)^^ t^-'[^~tf-'@{[K,}■-\-^^dt (11) 

(min{fR{a),5H((3)}>0; min{9^ (b) , <K (d)} > 0] 



Remark 



(I) For given sequence {K^ji^gj^^our definitions would reduce to known or new extensions of the Gamma and 
Beta functions. Recall that the asymptotic behavior of Kummer's confluent hypergeometric function 
at infinity has the form 



(D (a;c;z) = iF, (a;c;z) = ^e^z"-'^ 

r(a) 



1+0 



$n(z) 00. 



It is clear that O {a;c;z) is a special case of 0(ki,;z) with sequence | -r-r" l • And our Gamma 
function V^^^ and Beta function S'-'^^' are redi^ced to (6) and (7) which has been studied in [2j]. We 



can also replace iFi (a;c;z) with exp(z) to get more basic extension of Gamma and Beta function, 
(see, for details, [7|) 

(II) Function sj,^^' (a, |3) is a further extension of Beta function sj,'^'' (a, (3) by introducing complex pa- 
rameter d (?l(d] > 0). This improvement may be essential since —j — -j-^ — — (^ ^^[f^i)^^ ^ ■ Note 

that if we set d = b we will get function (a, (3). In general, form —j — -f4t gi'^G us more 

flexibility for making transformations. 





By using bI^^^ (a, (3) we can get the corresponding further extension of Gauss hypergeometric function 
as follows. 

Definition 1.2. fl^J The Extended Gauss hypergeometric function 2^1'^'' is defined by 



2h 



oc^ , (X2 



13 



;z;b, d 



n=0 



(a2, |3i - az] n! 



(12) 



(|z| < 1;fR{c) >m{h] > 0;min{9^(b), 91(d)} > 0) 
If O (ki; z) — exp (z), we write 

Bb,d (0C2 +n, Pi - CC2) z" 



2r^ 



ai, a2 



13 



;z;b,d 



00 



Tl = 



Its integral representation can be easily derived. 



B (a2, (3i - (X2) n! 



Theorem 1.1. For the extended Gauss hypergeometric function 2^\'^^ '^ (<^1 > ct2; |3i ;z; b, d), we have the fol- 
lowing integral represetation: 



21-1 



:(Ki) 



(3i 



;z;b,d 



1 



1 -t 



1 -zt)-«' e Ki;--- 



B (a2, |3i - (X2) Jo 
<K{|3i) >m[a2) >0; «K(b),<n(d) > 0;b = d = 0, |arg (1 -z) <7t| 
To prove this integral representation, we first expand (1 — zt) as its Taylor series 

(zt)- 

n=0 



t 1 - t 



dt, 
(13) 



(1-zt)-"' 



It is clear that this series is absolutely and uniformly convergent. Thus, we can interchange the order of 
summation and than integrate out t to get the result. 



As outlined above, we have introduced the general idea of the development path of the subject, the 
extended Beta function and its application in generalizing some hypergeometric series. As we all know, 
general hypergeometric functions have important applications in mathematical, physical and chemical areas. 
Thus, quite naturally, we should not neglect the question that whether we could extend the definition of 
generalized hypergeometric function by using extended beta function (a, |3). Or whether we could 

prove that general hypergeometric functions are specific cases of the functions we get in the paper. The 
significance of such questions lies in the fact that a large number of practical problems require more complex 
hypergeometric functions of complicated coefficients. Fortunately, the answers are almost certainly "yes", 
inspired by the way how extended Appell's hypergeometric functions F2'^'' (a, |3i, |32;Yi )Y2;'<^)'y;b, d) are 
constructed in [Tsj . (2l| . In fact, the method can also be applied to ratiocinate a new kind of multivariate 
hypergeometric function, i.e., F^ j^^^ j (a, (3i , • • • , (3r;yi > • • •Tt;''i > • • • )Xr;b, d), with its large numbers of 
formulas, which will be analyzed in section 3. 

Let us review some basic conclusions concerning with general hypergeometric functions. The generalized 
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hypergeometric function with p numerator and q denominator parameters is defined by [13, p. 27] 



Fq (oci,--- ,ap;|3i,--- ,|3q;z) = pFq 



V q 



ai 
Pi 



a. 



p . 



k=0 



(14) 



r(|3i)---r((3q) ^ r(ai +k)-..r(ap +k) z'< 



r(ai)---r(ap) f- r 



k=0 



:i3i +k)---r((3q +k) k! 

(15) 



(ai,|3j eC,|3j ^0,-1,-2,--- ,1=1,... ,p;j = l,... ,q) 

which is absolutely convergent for all values of z G C if p < q . When p — q + 1 , the series is absolutely 
convergent for |z| < 1 and for |z| — 1 when (^^jLi Pj — ^1=1 > 0, while it is conditionally convergent 

for |z| = 1 (z ^ 1 ) if — 1 < $H (^jLi Pj ^ ^1=1 '^ij — 0- More detailed information may be found in [2,|3]and 
2I. 



In what follows, we consider the extended generalized hypergeometric function pFq'^'' (ai , . ■ . ap; (3i . . . |3q;z;b, d) 
defined by 



pTq 



Pi 



Pq '^'^''^ 



^ ^b!d (o^i+i + TTt, Pi - gj+i ) z 



Tn.=0 



j = l 



B (aj+i, Pj - Kj 



+1 ^ 



°° ^ <^'(aj+m,pj^aj) z"- 
m=oj=i B(aj,Pj-aj) m!' 



m=0 



1 



(Pi 



1 /r <^'(ai+m,p,+j-aj) z"^ 
(Pr)mfJ B(aj,P,+j -aj) m!' 



|z|< 1; p = q + 1 
p = q 

r = q-p, p< q. 
(16) 



If we set b = d = 0, then 



F'^i' (ai,...ap;Pi . • . Pq;z;b, d) = pFq (oci , • . • , ap; Pi , . • . ,pq;z). 



P' q 



If p = 2, q = 1 then pFq*^^' reduces to extended Gauss hypergeometric function (12). 



The paper is organized as follows. 

In section 2, we will prove Euler type integral representation and expound some properties of new 
extended generalized hypergeometric function pFq'^^'. By using Ramanujan's Master Theorem, we can get the 
Mellin-Barnes type integral representation for extended generalized hypergeometric function pFq'^' ' . Then, as 
a special case of pFq'^'', more properties of extended Gauss hypergeometric functions such as differentiation 
formulas and recurrence formulas will be discussed in detail. In the last part of this section, we give an 
improved definition, definition 2.2. By doing this, we are not only able to prove a summation theorem 
deserving special attention, but also to establish the relations between modified functions and fractional 
operators. 

In section 3, we will give more integral identities and reduction formulas for functions such as extended 
Appell's hypergeometric functions F^'^'' (a, (3, (3';Y;x,'y; b, d), Fj''^' (a, P, (3';y, Y';x,'y;b, d) and Laucrella's 
hypergeometric functions Fj-, j^^^j, defined in |18||. In addition, we give a definition of a new type of multi- 
variate hypergeometric function, i.e.. 



^a,\ki) ("^fPi)--- >PT;Yi,---yr;xi,-.- ,Xr;b,d) = 



[\x^\ + ■■■ + \xr\ < 1; min{fR(b),fH(d)}> 0] 

and draw conclusions about it on the basis of the definition. Viewed as a a generalization of Ramanujan's 

Master Theorem and with its wide application in evaluating multidimensional definite integrals and the 
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Fcynman integral, the Method of Bracket is also used to prove the Mellin-Barnes integral representations of 
these multivariate hypergeometric functions. 

In section 4, we reestablish a Hilbert-Hardy inequality by using extended Gauss hypergeometric func- 
tions. It is remarkable that the extended Gauss hypergeometric functions can also be used to generalize 
other inequalities. 



2. Extended Generalized Hypergeometric Functions 

We first give a formal definition and some remarks for extended generalized hypergeometric functions. 

Definition 2.1. For suitably constrained (real or complex) parameters aj, j = 1, • • -p; (3i, i = 1, • • ■ , q, we 
define the extended generalized hypergeometric functions by 



■■■ "^"P -z-b d 

(3l ••• (3q '^'^''^ 



(18) 

Remark 



(I) If the numerator parameter ai in the first series of (18) is zero or a negative integer, the series 
terminates, i.e., 



q+1 Tq 



|3l •■• |3q '^'■'''^ 



^0 M B(aHi,P)-cxHi) m!- 



And we should make sure that no denominator parameter |3i , • • • , |3r is allowed to be zero or a negative 
integer in the third expression of (18). 

(II) The convergence of series pFq'^' ' with different parameter p and q can be obtained by using Comparison 
Test and the convergence of known corresponding generalized hypergeometric functions pFq . Their 
own convergence region is therefore clear. 

(III) Note that for usual generalized hypergeometric function p Fq ( ai , • • • , ap ; (3 1 , • • • , (3 q ; z) , its numerator- 
parameters ai , • • • ap and denominator-parameters pi , • • • , |3q can be arbitrarily rearranged without 
any effects on its value. But under the setting of our definition, we need to point out that the expres- 
sion shown on the right-hand side of (18) is closely associated with the order of the parameters given 
in the square bracket on the the left-hand side of (18). 

(IV) Although the general definition seems complicated, it will be very convenient in practice. 



We will first discuss in detail some fundamental properties of function (18) and then in the final part of 
this section a modified definition will be presented. 



2.1. Some basic results about the extended generalized hypergeometric function 

What we most concerned about is the Euler type integral representation of our new function (18). The 
following theorem aims to demonstrate that the form of the Euler type integral representation of pFq'^'' is 
very similar to that of the Euler type integral representation of pFq . 
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Theorem 2.1. //p < q + 1 , (|3q ) >5H(ap] >0 andm{h) > 0, m{d) >0, we have 



r(Ki) 



|3i 



;z;b, d 

Pq 



r(ap) r(|3q - apj Jo 

p(Kl) 



■■■ '^i'-^ -zf b d 



nl ^ 'I 

Ki;-- - 



t 1 -t 



dt. (19) 



Eq (19) also holds when b = d = 0, proided that jarg (1 — z) < 7t|. 
Proof 

We need to verify that formula (19) holds for three different expressions of pFq'^'' (ai , • • • , ap; (3i , • • • , pq;z; b, 
given in (18), respectively. Consider the case p = q + 1. Note that 



'^b'^d ("^q+1 + Tg, |3q - gq+i ) ^ F ((3q) 

B(0Cq + l,|3q -aq + l) r(aq + l)r(|3q - ttq + l ) Jq 



t 1 -t 



Substituting this in expression (18), we get 



q + 1 Tq 



ai 

Pi 



O^q + l 
|3q 



;z;b,d 



m=0 



B(aj+i,(3j -aj+i) r(aq+i)r(|3q -aq+i) 



dt 



r(aq + l) r(|3q - aq + l) Jq 



t 1 -t 

e( Ku 



m! 

b d 



t 1 -t 



X 

m=0 

r(pq) 



B(aj+i,(3j -aj+i) m! 



r(0Cq + l) r(|3q -aq + l) 



Jo \ t 1 — t 



X aF 



Kl) 



q ' q-l 



ai 
|3i 



p,'^'' ;zt;b,d 

Pq-l 



dt. (20) 



For p = q , we have 



r(Ki) 
q "^q 



ai 
|3i 



r'' ;z;b,d 

Pq 



r(|3c 



r (aq) r (|3q — aq 



Pq-ttq-l 



X iF'"'' 
^ q-l "^q-l 



q J I W-Jq ^ '^q) Jo 
CC^ ■ 
Pi • 



t 1 -t 



^^-^ ;zt;b,d 

Pq-l 



dt. 



It is clear that this relation is valid for p < q. This completes the proof. 
Remark 

A multidimensional case of Euler type integral representation of (20) is given by: 



□ 



q + 1 Tq 



ai ••• aq + l 



;z;b,d 



0^ 



r(Pj] 



(aj+i)r(Pj -aj 



+1 . 



b 

Ki;-— 



tj 1-ti 



(1 - tit2---tqZ) dti •■•dtq, 



which follows from the repeated application of the functional equation (20). If we set b = d = 0, then above 
equation reduces to the form given by 12, p. 132, Eq.(4.2)] (} 

Theorem 2.2. Fory < q + 1, we have the following differentiation formula: 



dz- r ^ 



ai 
Pi 



p!; '^'^''i 



pTq 



(P 



6 



(P 



qJn 



ai +n 
Pi +n 



ap + n 

pj+n'^-'^''^ 



(21) 



Proof 



Taking the derivative of q+i Fq'^^'with respect to z, we obtain 



_d_ 
dz 



q + 1 Tq 



(3l 



_d_ 
dz 



TR = 



+ 1 , 



V 1 TT ^b.d +m,(3j -gj+i) z 



m! 



m-1 



Tn=l 



3 = 1 



B[aj+i,(3j -aj+i) (m-1) 



Replacing m — > m + 1 we get 



dz 



— < q + l Tq 



(3l ••• (3q '^'^''^ 



2_ i'^iJm+i 11^ 



m=0 j = l 

rq „, . _ oo 



m! 



ai + 1 
|3i +1 



(cxj+i + 1, (3j - aj+i) m! 



nj=i Pi m=o 



'^''+' +^ -z-b d 

(3q+1 '^'^''^ 



Recursive application of this procedure gives us the general form 



dz- 1 ^ 



gi • • • gq+i , , 

PI ••• Pq 



_ (ai)^---(aq+i)^ 



Similarly, we can prove this result for the case p < q . 
For p = 2 and q = 1 we get the following corollary. 



ai +n • • • aq+i +n 
|3i+n ••• |3q+n '^'''''^ 

□ 



Corollary 1. 



d- 
dz" 



2F 



ai,a2 



;z;b,d 



(Pi). 



(x^ + n, a2 + n 
Pi +n 



;z;b,d 



(22) 



Remark 

If we take Q (ki;z) — iFi (oc; P;z) ,b = d, then equation (21) reduces to 



d- 
dz^ 



(a,P) 



a, b 



1-^1 b 



a + n, b + n 
c + n 



;z 



given in [22, Theorem 3.3]. <) 

In order to derive the Mellin-Barnes type contour integral representation of pFq'^'' (ai , • ■ • , ap; pi , • • • , Pq;z;b, d) 
we need to introduce the following well-known theorem. It is widely used to evaluate definite integrals and 
infinite series. 



Theorem 2.3 (Ramanujan's Master Theorem [l||). Assume f admits an expansion of the form: 



k=0 



k! 



Then, the Mellin transform of f is given by 



x'-V(x)dx = r(s)A(-s). 



By means of Ramanujan's master theorem we obtain: 
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Theorem 2.4. We have the following Mellin-Barnes type integral representation of function (18), namely, 



'^^ ■■■ "^"P -z-b d 

^ B[^^^(aj+1 -s,(3j-aj+i)r(s)r(ai -s) 



i n 



r(ai) 



(s) (-z)-^ ds, 



"z] ds, 



_L / A <d («j - s, Pj+r - cxj) r (|3i ) ■ ■ ■ r (p,) r (s] 

ZmiL^i B(aj,Pj+,-aj) T (Pi - s) • • • T (p, - s] 



p = q + 1 



p = q 



(23) 



-z) ' ds, r = q - p, p < q. 



where L is a Barnes path of integration, that is, L starts at —ioo and runs to +ioo in the s-plane, curving if 
necessary to put the poles ofr[a-\ — s) to the left of the path and to put the poles o/ r(s) to the right of the 

path. 

Proof 

The result follows directly by using Ramanujan's Master Theorem and inverse Mellin transform. □ 

In section 3, a Multidimensional Ramanjuan's Master Theorem will be introduced to obtain double and 
multiple contour integral representations for extended Appell's and Lauricella's hypergeometric functions. 

2.2. More properties about Gauss hypergeometric function 

Although many different kinds of extensions of Gauss hypergeometric function have been given, we still 
know little about their sorts of properties. Therefore, we list some formulas and conclusions about extended 
Gauss hypergeometric function,of which some are originated from classical theories and others are drawn on 
the basis of recent researches. 

Theorem 2.5. The following transformations hold true for extended Gauss hypergeometric function: 
1. PfEtfF transformation 



2'-^ 



P 



;z;b,d 



ai , P 1 — ai . z 

<X2 ' 1 - Z 



;d,l) 



, (|arg(1-z)|<7t). (24) 



2. Euler transformation When 0(kv;z) =exp(z], we have: 



ai, a2 
Pi 



;z;b,d 



^g-(l-z)b-zd(l _, 



iri 



Pi-ai P1-CX2 b ■ 
Pi 1 — z 



(largd -z)|<7t, $K(z)<l] 



(25) 



Proof 

1. The proof of transformation (24) is direct. Since we have [1 — z(l — t]]^" — (1 — z)^°' (l + jz^t) 
then, by replacing t ^ 1 — t, we obtain 



(Kl) 



' ;z;b,d 
Pi 



1 



B {(X2, pi - OC2) 

{^-z)-^ 



B [oil, Pi - Oil] Jo 



1^ (1 _t)«-i tP'-^-i [1 _t)]-«^ e (Ki;-y^ - ^) 



1 -z 



dt 
dt 



= 0-z) 



-a t:(ki) 
2^1 



ai , p 1 — ai . z 



;;d,b 



oci '^-z' 

2. When Q (ki,;z) = exp (z), from Pfaff transformation (24) we have 



ifi 



a^,al 



P 



;z;b, d 



(1-z)-"'2Fi 



ai,Pi -a2 . z 



Oil 



' 1 -z 



;d,b 



The rest of the proof will be presented in the section 4 since we need more extra results built in that section. 
□ 



We also give a somewhat complicated differentiation formula for extended Gauss hypergeometric func- 
tion. 

Theorem 2.6. One has 



+TI-1 t:(<i) 
2'"l 



<x^,o^2 



(3 



;z;b,d 



(3 



;z;b, d 



(26) 



Proof 

This formula can be obtained by direct computation, 



_d 
dz 



ai,a2 



;z;b,d 



(^2+n,(3i -a2)z^+"' 
B [(X2y (3i - ai) n! 



.n=0 



2_ (cxi)^ (ai +n) 



n=0 



B (a2, (3i - Oil) 



n! 



(Note that (a + n) (a)^ = a(a+ 1)^) = aiz"'"^ ^ (ai + 1)^ 



n=0 



B (a2, Pi - a2) Ti! 



aiz«i-'2F 



1 irtKi) 



oci + 1 , a2 
Pi 



;z;b,d 



dz2 



;z;b,d 



= Y_ (o^On ('^i +n+ 1) (oci +n; 



S^(a2+n,|3, -a2)z-+«'-i 



n=0 



ai (a, +l)z"i2F 



(Kl) 



ai + 2, 0C2 
Pi 



B (a2, Pi - 0C2) nl 



;z;b,d 



In the procedure, we have used the fact that (a + n + 1 ) (a + n) (a)^ = a (a + 1 ) (a + 2)^. Recursive appli- 
cation of this procedure give us the general form (26). □ 
A class of recurrence relations for Gauss hypergeometric functions 



2F1 



ai ± n, az 
Pi 



, 2F1 



ai , a2 ± TL 
Pi 



and 2F1 



ai,a2 
Pi ±n " 



(neNo) 



have been established in (20| as a useful tool to find some new recursion formulas for Appell hypergeometric 
functions 

F2 (cT, ai, a2 in; |3i, P2;x,'y) and F2 (cr, ai, a2; Pi, P2 ±n;x,y) . 

It is natural to consider whether these results are true for our new definitions. The following theorem 
gives us a positive answer. 



Theorem 2.7. The following recurrence relations hold true for the extended Gauss hypergeometric functions 

2F1 . 

1. Recurrence relations for 2F1'^'' (ai ± n, 0C2; Pi ;z;b, d) 



2F 



(Kl) 



ai + n, 0C2 



P 



2F 



(Kl) 



ai — n, a2 



P 



;z;b, d 



;z;b,d 



= 2F 



(Kl) 



0Ci,0C2 



P 



2F 



(Kl) 



ai,a2 



P 



;z;b,d 



;z;b,d 



0C2Z ^ 



(Kl) 
1 



lc=l 



ai + n — k + 1 , a2 + 1 
Pi +1 



U2^ V- 
' k=l 



(Kl) 



oci — k -h 1 , a2 + 1 
Pi +1 



|zl < 1; neNo) 



;z;b,d 
(27) 
;z;b, d . 
(28) 



2. Recurrence relation for 2F]'^'' (ai,a2;Pi +n;z;b,d] 



2F 



(Kl) 



ai,a2 
Pi +n 



;z;b,d 



(Pi 



(Pi -^2), 



k=o 9 



k / n \ (a2) 



k t:(ki) 



^ J (Pi) 



^2F 



ai , a2 + k 
Pi +k 



;z; b, d 



• (29) 



\z\ < 1; neNo) 



3. Recurrence relation for 2^1'^'^ (cci , ai + n; |3i ;z;b, d) 



(Kl) 



' ;z;b,d 
pi 



((3i -a2)2n 

(Pi -a2)i^(a2)n 



(Pi 



ai , a2 + i + n , , 



:iz| < 1; neNo) 



(30) 



Proof 

1 . By means of Euler integral representation of the extended Gauss hypergeometric functions given by 
(13), we have, by writing (1 -zt)""'""" = (1 -zt)""^^''-^ -zt(1 -zt)""^"''"\ 



2^ 



(Kl) 



a^ + n, cx.2 
Pi 



;z;b,d 



r(Pi) 



r(a2)r(|3i -az] 

zr(|3i) 

r(a2)r(|3i -a2) 



|\oc2-^ (i_t)Pi-«2-l (l_2t)-«i-n-l0j^^^._b__^^ 
^0 



dt 



=2Fr 



Replacing n by n — 1 , we find that 



_ p(Ki) f ai + n + l,a2 



Pi 



;z;b,d 



ai + ny dz 



b d 
"t " 

ai +n,a.2 



Pi 



;z;b,d 



2f; 



(Kl) 



ai + n, 1x2 
Pi 



;z;b,d 



= 2F 



(Kl) 



ai +n—^yOC2 



+ 



;z;b,d 



ai + n — 1 y dz 



2F 



(Kl) 



ai + n - 1 , a2 , , 
Pi 



(31) 



Applying this identity recursively, we obtain: 



2f; 



(Kl) 



ai + n, a2 , , 

Pi 



2F 



+ z > 

^ fYi -I- 1 



ai,a2 
Pi 
d 



;z;b, d 



k=1 



= 2F'^^^ 



ai + n — k dz 



2f; 



(Kl) 



ai + n — k, a2 



Pi 



;z;b, d 



ai,a2 



P 



;z;b, d 



lc=l 



ai + n — k + 1 , a2 + 1 
Pi +1 



;z;b,d 



(32) 



To prove the second formula, we still begin with the integral representation of extended Gauss hypergeometric 
function. Writing (1 -zt)"'"^"^^ = (1 -zt)-^'^'-''+^^ -zt(1 -21)"''^'-^+^^ we have 



2FI 



Kl) 



ai — n, a2 



Pi 



;z;b,d 



r(Pi 



r(a2)r(|3i -a2) 



r(a2)r(|3, -a2) 



f (1 __tj|3,-«2-l (1 _^t)-(«i-n+i)0/ _ d 

Jo \ t 1 — t 

{^-tf'-'^'-' (1 -zt)-'"^-^+'^ e (^Kv-^ ~ 



dt 
dt 



:(Ki) 



= 2F, 

Replacing n by n + 1 , we find that 



ai - n + 1 , a2 , , 
;z;b,d 

Pi 



ai — n y dz 



£ 2f; 



:(Ki} 



ai - n, a2 , , 
;z;b,d 

Pi 



(33) 



2f; 



(Kl) 



ai — n, a2 , , 
;z;b,d 
Pi 



2F 



(Kl) 



ai — n — 1 , a2 , , 
R ;z;b,d 
Pi 



— z 



dz [ n — ai + 1 
10 



1 p(Kl) 

2^1 



ai — n — 1 , a2 



;z;b, d 



(34) 



Repeating this recurrence relation n times and applying Corollary 1, we obtain: 



2f; 



ai - n, a2 , , 



(Kl) 



ai - In, 0C2 , , 



|3 



k=1 



ai — n — k + 1 , a2 + 1 
Pi +1 



;z;b,d 



(35) 



The identity (28) follows directly from (35) upon replacing ai by ai +Ti (n G No). 
2. Using Euler integral representation (13) and the fact that 



(1-t) 



3i+n— a2 — 1 



= (1 -t)P^-"'-^ (1-t)^ = (1-t) 



Pi-a2-l 



we find 



2F 



(Kl) 



ai,a2 
Pi +n 



;z;b,d 



r(pi +n) 



r(a2)r(Pi -az+n] 



k=0 



lc=0 



■.ai +lc— 1 



(i-t) 



Pi +lc-(a2+k)-l 



xn-zt)«-^e(K,;4- ^ 



r(Pi +n) 



r (a2) r (Pi — (X2 + n 



t 1 -t 

TL 



dt 



n \ , r(a2+lc)r(Pi -a2) ^(k,) 
I.~1J FTT; — TT^ 2ri 



(pi -a2). 



k=0 



k=0 

iPiJk 



r(Pi +k) 



ai , a2 + Tc 
Pi +k 



;z;b,d 



ai , a2 + k 



Pi +k 



;z;b, d 



This completes the proof of Eq (29). 

3. The demonstration of assertion 3 is based on the following elementary expansion 



i=0 



(36) 



i=0 



Substituting expansion (36) in the integral representation of extended Gauss hypergeometric function (13) 
we have 



2F 



(Kl) 



ai , a2 + n 
Pi 



;z;b, d 



r(Pi) 



r(a2 + n] r(Pi - a2 -n) Jo 
b d 



i-ai+n— 1 



(1 _t^Pl-«2-n-l _^^y 



X e Ki; 



t 1 -t 
r(Pi 



dt 



r(a2)(a2)^r(p 



-T t. (-^^i / t«^+-+-^ (1 - t] 

1 - a2 - n) ^ Vo 



Pi-ai-l 



X (1 -zt)-«' G(^KU- 

^ npo 

r(a2) (a2)^r(Pi -a2-n; 



t 1 -t 

n 



dt 



1=0 



r(a2+i + n)r(pi -a2) 
r(p, +i + n) 



X2F; 



(Kl) 



ai , a2 + n + i 



Pi +n + i 

r(pi -ai) 
(a2)„ r(Pi -a2-n) 



;z;b,d 



f^2Ji+n^p(Ki) 



i=0 



(P 



ai , a2 + n + 1 , , 
Pi +n + x 



And equation (30) follows by using the following property 

r(pi -a2] 



r(Pi-a2-n) - - + = (Pi-a2),' 



□ 

RemEirk 

For classical Gauss hypergeometric function, the parameters ai and ai are symmetric, i.e., 

2F1 (ai, aa; |3i;z) = 2^1 (ai, ar, (3i;z) . 

Thus we only need to establish the recurrence relations for one of them and the results can be valid for both 
ai and ai- However, the numerator parameters for extended Gauss hypergeometric functions possess no 
such symmetry. 

2.3. A further generalization and its connection with fractional calculus 

Our definition of extended generalized hypergeometric function can in fact be generalized to the following 
form: 

Definition 2.2. For suitably constrained (real or complex) parameters aj, j = I,-- - ,p; Pi, i = I,-- - ,q, 

we define the extended generalized hypergeometric functions by 



(ai,ki) 
|3i, 



P q ) 



m=0 



B(aj+i,Pj -aj+i; 
B (aj, |3j — aj) m! ' 



m!' 



m=0 j = l 



|z|<l; p = q + 1 



p = q 



(37) 



l,m=0 



(|3i 



where the new introduced parameters kj , j = 1 , • • ■ , p are non-negative integers. 

Obviously, (37) reduces to (18) whenever kj = 1, j = 1, • • • ,p. To illustrate its advantages we first 
consider the following function: 



(ai,1),(a2,k2) 
|3i 



;z;b,d 



B(a2,|3i-a2) n! 



n=0 



(38) 



Its integral representation can be written as 



(Kl) 



(ai,1),(a2,k2] , , 
|3i ' 



1 



B(a2,|3i -oti) 



t«^-i [^-tf (1 -zt''^)""^ (^Kv-j - dt. (39) 

9l(|3i) >fH(a2) >0;fH(b), 91(d) >0;b = d = 0,|arg(1 -z) <7r| 



12 



As an example, let us compute the case for ki = 2 and z—^, 
.1 



r((3 



r(a2) r((3i - az) Jo 

r(Pi) 



dt 



r(a2)r(|3, 




r{pi 




r(a2)rOi 


- ai) 


r(Pi 





/' t«^-i (1 -t)"^'-"^-^ (1 -t)-«' (1 +t)-"' e (ki;-- - 
Jo \ t I — t 

oo 



dt 



r(a2)r(|3, -az) 



k=0 



-ai 
k 



(1 -t) 



G Ki;-- 



t 1 -t 



dt 



(-1) ('^l)k r , 1 a 

■Bl \' (a2 + k, Pi -a2-a,. 



k! 



k=0 

r(Pi -a2-ai)r(a2) 



r(a2)r(|3, -a2) r((3i-ai: 



k=0 



gj,^^' (a2+k,(3i -a2-ai) ("1)'' 
^ B (a2, |3i - a2 - oci ) 



k! 



Pi - ai 



;-i;b,d 



r(pi -a2)r(p, -ai; 

We conclude our result as the following theorem: 

Theorem 2.8. If^[^^] >m[(X2] >0 and<n(Pi - a2 - ai ) > 0, then 



2F 



(ai, 1) , [(x.2,2] 
(3i 



;i;b,d 



r(Pi)r(Pi -(X2-(x{_ 
r(Pi -a2)r(pi -ai: 



2F 



OCi, 0C2 
Pi - OC] 



;-1;b,d 



(40) 



Theorem 2.8 is in fact a generalization of the following well-known result [i, p.117, Theorem 2.3]: 



3F2 



n b b + 1 
2 ' 2 .1 
c c+1 ' ' 
2' 2 



r(c)r(c-a-b) 

r(c-b)r(c-Q) 



2F1 



a, b 

c — Q 



-1 



(41) 



It also show some intrinsic relationships between our extended hypergeometric function and the ordinary 
generalized hypergeometric functions. 

In general, we have for definition (37) the following integral representation, which is a generalization of 
US, p. 104, Theorem 38]. 

Theorem 2.9. //p < q + 1, 5H(Pq) > 5H(ap) > 0, kj e Z+ U{0},j = I,--- ,p and none o/Pi,...,Pq is 

zero or a negative integer, we have ( inside the region of convergence ) 



t:(ki) 



(ai,ki; 
|3i 



;cz'^p;b, d 



r(Pq)zi-p^ 

r(ap)r(Pq -ap] 



X T,-iF 



(Kl) 



P-1 ""q-l 



X O K;; 



(ai,ki ) 
|3i 

bz dz 

t z-t 



Pq-1 



dt. 



(42) 



Proof 
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Denote the right-hand side of (42) by A (z). As indicated before, we just need to prove the case that 
p — q + ]. Let t = zv in A (z) . Then we have 



Afzl = 



r(|3q)zP^-i 



X nF 



(oCn, kn 



r(0Cp)r(|3q -OCp) Jo 

(Ki) r (ai,ki) •■ 

r((3q)zP^-1 y 

Jo 



''"""^ ;c(zv)''^+' ;b,d 



1 Ki; I dv 

V 1 — V , 



q-1 c(k 



n 



^b.d + TT^' Pi - ) (^^ 



-q +1 



i=i 



e K 



dv 



ra! 



V 1 — V 



"'"M B(aHi,|3i-aj+i) m! 



— Z q+1 Fq 



(ai,ki) ••• {aq,kq) (aq+l,kq+i) .^,^k,+,.^^^' 
Pi ••• |3q-l (3q ' ' ' 



□ 

Theorem 2.9 embody the significance and advantages of our modified definition (37) of extended gen- 
erahzed hypergeometric function. It not only enables us to rebuild a lot of classical results of hypergeometric 
function in a new way but also provide some new important interpretations. Here is one of them. In [l8| . 
H. M. Srivastava definite the following fractional derivative operator: 



V 



|a.,b,d 



\l (z - 1)-^-^ e (ki; - ^) f (t) dt (91 (^) < 0) 



(43) 



(m- 1 < 9^ (m-) < Ta(m e N)) 



The path of integration in (43) is a line in the complex t— plane from t = to t = z. 
By using operator (43) we can rewrite (42) as 



pTq 



(oci,ki; 
(3i 



r (Pq) ^1-|3„-p-(|3q-«p),b,d f.ap-l p(Ki) 



(ai,ki; 
Pi 



^^'v-^^^v-^^ ;ct'^p;b,d 

Pq-l 



(44) 



3. Extended Appell and Lauricella Hypergeometric Functions 

In this section, we focus on some extended multivariate hypergeometric functions, i.e., extended Appell's 
hypergeometric functions and extended Lauricella's hypergeometric functions. 

For the sake of clarity and easy readability, we may first study the properties of extended Appell's 
hypergeometric functions and then we can view the extended Lauricella's hypergeometric functions as a 
further generalization of Appell functions. 

Definition 3.1. The extended Appell's hypergeometric functions 



(a,(3i,(32;Yi;x,'y;b,d) and F^"'' (a, |3i, |32;Yi,Y2;x,y;b, d) 



in two variables are defined by 



F5^^'(a,(3i,(32;Yi;x,y;b,d)= ^ (|3i)^((32)^ 

m,n.— 



bI"'^ (a + m + n,Yi -cc) x"^y^ 



B (a,Yi - tx) 
max{M,|y|}< 1 ; min{$n (b) , (d)} > 0) 



m!n! 



(45) 
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.M, a a u ^^ \^ f ^ ^1,^^' (Pi + m,yi - p, ) b{,^^' (Pa + n, - Pi) x-y- 

^•^+^ B(Pi,Yi -Pi)B(P2,Y2-P2) m!n! 



F2 (a, Pi,P2;Yi,y2;x,y;b,d) = Y_ 

m,n=0 -u-.w. i-i. - V,--, 

(46) 

(|x| + |y|<1;min {91(b), fH(d)}>0) 

The following Euler type integral representations for Extended Appell hypergeometric functions F^'^^^ 
and F2''^^ are easily obtained and their proofs will be omitted. 

Theorem 3.1. For the extended Appell functions defined by (45) and (46), the following integral represen- 
tations hold true: 

f\^'^ (a,pi,p2;Yi;x,y;b,d) 

(min{9l(b),fH(d)}>0;max{|arg(1 -x)|,|arg(1 -y)|} < 7t;fH(Y) > fH(a) > 0) 
F^"'^ (a, Pi,P2;Yi,Y2;x,y;b,d) = ^ ^ 



B(Pi,Y, -Pi)B(P2,Y2-P2) 
/ / —r. Ki; Ki; q dtds. (48) 

yoyo (1-Xt-ys)" V t l-ty/ V S l-Sy/ ^ ^ 

(min{$K(b),in(d)}>0;|x| + |y| < 1;$K (yi ) > 9^ (Pi ) > 0,1H (Y2) > 91 (P2) > 0) 

By using the integral representations of F^'^^' (a, P, P'; Y;x,y;b, d) and F2*''' (a, P, P';Y,Y';x,y;b, d), 
we can derive some transformation formulas for them. 

Theorem 3.2. For the extended Appell hypergeometric function f\'^^\ we have the following transformation 
formula: 



F^^^' (a, p, p';Y;x,y;b, d) = (1 -x)"'' (1 -y)"'' F^^^' (^a, p, P';y; — j-, d,bj . (49) 



Proof 

Let us put t = 1 — s. Then 



F^^^' (a,p,p';y;x,y;b,d) = — -! . f s^""-' (1 -s)'^"' (1 -x(1 -s))"'^ (1 -y (1 -s))" 

B (a,Y- a) Jo 



xe(K,;-:A__£)ds 



B(a,Y-cx) 



/ s^-*-^ (1 -s)""^ (1 -x + xs)"'^ (1 -y +ys) 

yo 



/^s--O-s)-(l-^s)- ( 



1 2 s 

B(a,Y-a) Jo " ' V x-^"J \ y-1 

x0(Ku-^-y^) ds 



such that 



Ft^^) (a, p, p';Y;x,y;b, d) = (1 -x)-^ (1 -y)^' F^^^' (^y - a, p, P';y; d,b 

□ 

Remark 
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The classical theory of Appell hypergeometric functions have found that integral 

-1 

u"-i (1 -uf-'^-' (1 -ux)-^ (1 -uy)-^' du 



keeps unchanged under the following five transformations [24, p. 195] 

u = 1 — V, 

V 1 — V 



1-1J+V1J> l-vx> 

V 1 — V 



-vy 

However, only the first transformation u = 1 — v enable us to get a simple extension. ^ 



Theorem 3.3. Forh — d we have the following transformation formulas for Appell hypergeometric function 

pf^'' (a, (3, (3';Y,Y';x,y;b) = (1 -xJ-^'F^^^' (^a,y - (3, (3';y,y'; ^;b) , (50) 

F^^'' (a, (3, (3';y,y';x,y;b] = (1 -yj-^F^^'' (^a, (3,y' - (3';y,y'; , (51) 

f'^'' (a, 13, |3';y,y';x,y;b) = (1 - x - F^^^' ( a,y - |3,y' - (3';y,y'; -— -— L;b) . (52) 
Proof 

By using the double integral representation of F2'^'' with b = d and the transformations 

(a) t = 1-t', s = s'; 

(b) t = t' s = s-1; 

(c) t = 1-t', s = 1-s'; 

we can directly obtain the transformation formulas for F2'^'' (a, (3, (3';y,y';x,'y;b), respectively. 
Note that formula (52) is also valid for more general case, namely, 

F^^^'(a,|3,|3';y,y';x,y;b,d) = (1 -x-y)-«F^^'' ('a,y - p,y' - (3';y,y'; -— -— L; d, b ) . 



(53) 
□ 

3.1. Other Integral Representations of Appell' s Hypergeometric Functions F^^'' andY^-^^^ 

We now start to establish the Mellin-Barnes type integral for Appell functions F^^'' and F2^^'. The 
tool we apply here is called Method of Bracket [1], which can be view as a multidimensional extension 
to Ramanujan's Master Theorem. This method enable us to evaluate multidimensional Mellin transform 
directly and efhciently. In this part, we will prove the following theorem with this method. Here we will put 
focus on the proving instead of explaining the method in detail. Because the method is easy-understanding, 
we will give a brief introduction of it when proving Theorem 3.4. More detailed information may be found 
in Q and 

Theorem 3.4. For suitable parameters, the Mellin-Barnes integral representation o/F^'^'' is 
F^^^'(a,(3i,(32;yi;x,y;b,d)^ .J.irm \^ri\ri 1 P^^b^d -S2,yi -a) 

(2711)-^ r(|3i) r(|32) r(oc) r{yi - a) 

X r (|3i - SI ) r ((32 - S2) r (si ) r (S2) (-x)-'^' (-y)-^^ dsi ds2. (54) 

Proof 
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We first take the Mellin transform of (ct, (3i , PilTi ; — x, —y,h, d), i.e., 

/ / x--iy--^Fi-'(a,Pi,P2;Ti;-x,-y;b,d)dxdy= ^ (Pi U (Pa)^ ""^"^ (a + m + n,yi-aj 
Jo Jo m Tt-0 di^a,yi— aj 

y x^'+'^-^dxJ ^ 4)m4)nf(m,n)(si +m)(s2+n) (55) 



-1]^ f-1 



r(m+l) r(n+1 
where 



m,n=0 



f(m,n) = ((3i)^((32) 



^b!d tcx + m + n,yi -g) ^ r(Pi +m)r(p2+n) g^^^^ (a + m + n,Yi - «) 



B(a,yi-a) r(|3i) r(|32) B(a,yi-a) 

'pTTi = -pT — r~rT>'PT 



r(m + i)'^'' r(n + i)' 

/'CXD /'OO 

Jo 

The symbol ^m. and 4)n are called the indicator of m and n, respectively fll'. p. 53, Definition 2.2]. Then 
use the method of brackets Rule III [1] we have 

x^i-iij^^-iF5^^'(a,(3i,(32;yi;-x,-y;b,d)dxdy = ^^^r(-m*)r(-n*)f(m*,n*) 



r((3i -si) r(|32 - S2) d ('^^ si - S2,yi - a 



-r(si)r(s2], (56) 



r(|30 r(P2) B(a,yi-a) 

where (m*,n*) is obtained by vanishing the expressions in the brackets (si +Ta), {s2+Ti). In general, 
(Ta*,n*) is the unique solution to the linear system obtained by vanishing the expressions in the brackets. 
And the determinant of its coefficient matrix is denoted by det(A). Here, we can see (a) is only a symbol 
associated with divergent integrals Jq°°x'^^Mx. But by introducing such a symbol, we in fact reduce the 
evaluation of definite integral to solving a linear system of equations, which is the basic idea of the method 
of brackets. 

Now we can take the inverse Mellin transform of (56) to get 

(27ti)^ J-ioo J-ioo r(|3i) r(P2) B(a,yi-oc) 
X r(si)r(s2)x-'^'y-'Msids2. 

Theorem 3.4 follows from replacing — x and — y with x and y. □ 
The Mellin-barnes double contour integral representation for extended Appell's hypergeometric function 
F2'^'"' (a, (3i , (32;yi ,y2;x,'y;b, d) is given in the following theorem. We prove this just the same way when 
proving Theorem 3.4. 

Theorem 3.5. For suitable parameters, the Mellin-Barnes integral representation of ''^ 

F2 a, Pi,P2;yi,Y2;x,y;b, d = — — ^^FTH^TFT a \ria irr a irr 1 

ilnif r (|3i ) r (yi - |3i ) T [^2) T (y2 - ^2] V [a] 

/ f?{,^^'((3i -si,yi -p,)B(,^^'((32-S2,y2-(32)r(a-si - S2) T (si ) T (52) (-x)"^' (-y)-^Msi ds2. 

-ioo J —loo 

3.2. Recursion formulas for Appell hypergeometric function F^^'''' 

Recursion formulas for classical Appell hypergeometric function F2 have been studied in [2^ by using 
contiguous function relations of the classical Gauss hypergeometric series 2F1 . These recursion formulas can 
be used to evaluate the radiation field integrals. Now we try to establish the similar formulas for our new 
Appell hypergeometric function F2'^''. 
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Theorem 3.6. The extended Appell hypergeometric function satisfies the following identities: 

F^"'' (a,|3i,|32 +n;yi,y2;x,'y;b,d) 

' (,rjJ:',r,l„ i'-"'.fe''i'''(''.g-.fe+" + ';T,,T.+i + n;x.,;b,d, (57) 

(|x| + \y\< l;inm{9l (b) , (d)} > 0; n G No, a, |3i , Pi € C; y, , yi G QZ^ ) 

and 

F^*"'' (a,|3i,|32;yi,y2 + n;x,y;b,d) 

= 7;r^^^^Z(-1)''f k ) 7^F^'(«,Pi,(32+k;yi,y2+k;x,y;b,d). (58) 

(|x| + lyl < 1;min{9l (b) ,fH (d)} > 0; n G No, a, |3i , |32 G C;yi ,y2 G C\Zo ) 

Proof 

1. Using the integral representation (13) the following double integral representation 
F^*"'' (a, |3i,|32;yi,y2;x,y;b,d) ^ ^ 



B((3i,y, -|3i)B(|32,y2-P2) 



^0 



Ki;---- e Ki; dtds 



(l-xt-ys)'" V t 1-t; V s 1-s 

can be reduced to a single integral involving the extended Gauss hypergeometric function 2l^i'^^\ i-e., 

F^^^) (a,|3i,|32;yi,y2;x,y;b,d) = ——J — t^^-' (1 -t)^'-^'-^ (1 -xt)-"' 

t> lpi,Ti - PI J Jo 

The assertion (57) then follows by substituting (30) into integral representation (59) with (32 = (32 +tv. Thus 

F^"' ^ (a, p 1 , P2 + n; y 1 , y2; X, y ; b, d) 

■ tP'-i (1 -tr'-P'-i (1 -xt)-«2Fi^^' ( <^^^^+'' ;:r-^;b,d^ 
Jo ' \ yi ^ -xt J 



1 



B(Pi,Yi -P 

x0(|Kv;-^-y4^^ dt (60) 

1 (Y2-P2]2n MP 



B(Pi,Yi -Pi) (Y2-P2)^(P2 



nUJ2J^— (Y2)i+^'-:Jo 

X 2Fi-' f ^1+1+^ ;:r-^;b,d') e ("k^;-^ - dt (6I) 
= r ^^Vf'!l"^ ^(-n)^ /'^"^^+" F^-^'(cx,Pi,P2 + n + i;yi,y2+i + n;x,y,b,d). (62) 

(Y2 - p2)n lP2 (Y2]i+n X! 
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2. The proof of the second assertion of Theorem 3.6 foUows directly by substituting from the the 
assertion (29) of Theorem 2.7 into the integral representation (59) with 72 —yz+n, 



(a, Pi, P2;yi,y2 +'n.;x,'y;b,d) 



1 



72 + Tl ' 1 - xt 



;b,d 



X B Ki; 

' t 1 -t 



dt 



(Y2 



(yi-Pa) 



k=0 



lYr-Pi-l 



1^ y [yik Jo B(|3,,yi -|3i 



[1 -xt) 



X 2F 



2h 

(yi) 



Ki) / a, (32 + lc , y 



y2 + lc '1-xt' 



;b,d e Ki; 



1 -t 



dt 



[yi-M 



k=0 



(63) 



□ 

Remark 

The second assertion (58) of Theorem 3.6 is clear a extension of the formula given in [ifll, p. 549, 
Theorem 3] and the first formula (57) seem to be new. In addition, the results obtained here are also 
applicable to the extended Appell hypergeometric functions 

ri"'' (a)Pi)p2;yi +n,y2;x,'y;b,d) and F^"'' (a,Pi + n, (32;yi ,y2;x,'y;b, d) . 
It seems to be a little difficult to establish some recursion relations for the functions: 



F2'"' (a,Pi,p2;yi -n,y2;x,'y;b,d) and (a,Pi - n, (32;yi ,y2;x,'y;b, d) . 



3.3. Finite sum representation of the extended Appell's hypergeometric function ¥\^^ 

Recently, Cuyt et al. @ obtained a finite algebraic sum representation of the Appell hypergeometric 
function Fi in the case when the parameters a, (3i , (32 and yiare positive integers with yi > a. Motivated 
by their works, we show that their results can be extended to obtain a finite sum representation of extended 
Appell's hypergeometric function F,'^''. 

Before we state our results, we need the following important lemma. 

Lemma 1. 0/; Let a — -g^, P — -^-^ and X = (1 — ux)^^ , Y = (1 — uy)^V Then for any positive integers 
s and t, 



t-l 

j=0 



i + s - 1 

s-1 



s-1 



k=0 ^ 



1 



k vs— k 



(64) 



Theorem 3.7. For any non-negative integers s,t and x 7^ y, we have for |x| < 1, lyl < 1 

t-l 



Ft^^'(1,s + 1,t+1;2;x,y;b,d)=y^+i^ 



j=0 
s-1 



] + S 

s 



-1 1 



(t-j + i),i 



;y;b,d 



k=0 
t + S 



k + t 



.k+t+1 ^T^l 



(ki) 



(s-k+1),1 



;x; b, d 



-1)^X^1, 



t+s+l 



y2F 



(ki) 



2 ;y;b,d 



X2F 



(Kl) 



1,1 



;x;b, d 



(65) 



Proof 
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Put a—], Pi = s + l, |32 = t+ l and Y2 = 2 in the integral representation (47), we obtain 
fJ^'' (1,s + 1,t+ 1;2;x,y;b,d) = (1 -xu)"'''+^' (1 " ^^r"^^ ' © (^ki; - du. (66) 



By using lemma 1 we have 



-(t+i) s+i \- n + s 



j=0 



-(t+l-j) , ot+l 



— UX 



,-(s+l-k) 



Substituting (67) into (66) one get 
Fl^^'(1,s + 1,t+1;2;x,y;b,d) = a^+i^ 



(67) 



i=o 



1 ^ n + s 

s 



(3' / (1 -uy]"'*+^^'e I Ki;-- 



k=0 



t 

Z 

j=0 

k=0 
t-1 



u 1 — u 

(1 -ux) ' ' e K:; 

u 1 — u 



du 



du 



(t-j + 1),1 



;y;b,d 



1 



(s-k+l),l 



;x; b, d 



t-i 



j=0 
s-l 



j+s+1 



;y;b,d 



k=0 



k+t+l 2'"l 



(s-k + l),l 



;x; b, d 



t + s\ (-l)*xty'^ 



t+s+i 1 y^' 1 



2 ;y;b>d 



2 ;x;b,d 



□ 



where we have apphed the integral representation of extended Gauss hypergeometric function. 
Remark 

If we set b = d = , then (65) is reduced to the following finite sum representation p. 215, Theorem 
2.1a] 



Fi (1,s + 1,t+1;2;x,y) 



(y -X 
t-i 

-z 

j=0 



x) -y'* / s + t 



s+t+l 



1 + S 



-x)'y^ 



[ln(1 -x) -ln(1 ~y]] 

i-(i-y)'"'T 



(y-x)'^+'+^ (t-j) 



k=0 



1 - (1 -X 



.k-s 



[x-yr^^^ (s-k) 



(68) 



<> 



Remark 

This method can be used to find more complicated finite sum representations of extended Appell's 
hypergeometric series F^'^^'. For instance, other finite sum representations given by [8|, p. 215, Theorem 2.1 



b,c] can indeed be considered analogously in a simple and straightforward manner. 







It marks the end of the discussion of extended Appell's hypergeometric functions. Next, we will put focus 
on the Lauricell's hypergeometric functions. However, not all four Lauricella's functions can be generalized 
to their new forms. here, we will give a known generalized result. 
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Definition 3.2. fl^] Lauricella hypergeometric function fI^'j^^^j is defined by 



D,(Ki) 



(a, Pi,-- - ,(3r;Y;xi,--- ,Xr;b, d) 



(br) 



'^b'^d («: + mi-| hTar,y-" Qc) x 



mi! mr! 



ZTT- (69) 



TTi] , ■ ■ ■ ,m.r —0 

(nmx{|xi|,.-. ,|x,|}< 1;min{lR(b),lR(d)}>0) 

(rl 

Its Euler type integral representation of Lauricella hypergeometric function j^^j is stated in the 
following theorem. 

Theorem 3.8. For the extended Lauricella hypergeometric function Fj^'^^^j defined by (69), the following 
integral integral representation hold true: 

^d!(ki) ,PT;r,xi,--- ,Xr;b,d] 



r(y] 



r [a] r (y — a) 



t«-' (1 -tr-«-^ n n -^jt)-p' e (^Ki;-^ - ^) dt. (70) 



(9^(b),<H(d) >0;b = d = 0,max{|arg(1 - xi )| , • • • |arg (1 - Xr]|} < 7t; 91 (y) >fR(a) > 0) 
As a direct consequence of above theorem, we have the following summation formula 



D,(Ki) 



(a, (3i,--- ,Pr;y; 1,--- ,l;b, d] 



If we put Xi = X2 = • • • = Xr = X, 



r[(x)r(y — a) ' 



(a, Pi, - • • , (3r;y;x, • • • ,x;b, d] = 2^1"'' (a; Pi + ■ • • + Pr;y;x;b,d) . 



Pr) 



Theorem 3.9. One has 



xrJ:',.,.fa,-Ai,...,-A,;a+P;-i^-^,...,- 



(b - q) fr p q 



D-(-i) V"^' ' '■"'^ ' afi+gi' ' af, + g, 'b-a'b-Q 

where a, b G M (a < b), fi,gi,Ai e C, i = 1, • • • ,k,fR(b) > 0,$n(d) > 0, m{a) > 0, 9^(P) > and 

< 1. 



max 



(b-a)fi 



cifi + gi 



(b-a)fr 



afr + gr 



Proof 



Our proof mainly depends on the following expansion (l7| 
(f,t + g,) >=(afj+gi) ^ 

mi=0 J ^ J 3J 



(71) 



(|(t-a)fj|<|Qfj + gj|; t e [a, b]] . (72) 
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Opening up ri[=i (fjt + 9))"^' by using (72) and then integrating out t we have 
£ (t - a)'^-^ (b - if-' n (f^t + (-^•' - 

)=1 mi ,---mr=0 ^ ' ^ ' 

X y''(t-a]™+-+^'-Vb-t]P-^e(^Ki;-;^-:^^dt. (73) 
The last integral in (73) can be evaluate by applying substitution t = (b — a) u + a, namely, 

p q 



t-a b-t 

^a+P+miH l-mr-1 / a+miH hrrir-l f i _ , , A p-1 



dt 



= (b-a)«+P^"^^+-+^'-' / u<^+"^^+-+'^^-'(l-u)P-'G Ku-7 r du 

Jo V (b-a]u (b-a)(l-u)y 

= (b - a]«+P+'-'+ -+"^'-i b'-^ , (a + mi + • • • + lur, |3] . (74) 

1j — a ' b — a 

Thus, 

£ (t - a) «-^b - 1) n (f, t + gj e (k,; - ^) dt 

r oo 
j=1 mi ,---mT=0 

1 / (b-a)fi\'^^ 1 / (b-a)fr 



mi! V afi +91/ TUr! V afr + 9r 

r 

(b-a)«+P-^n('^fj + 9j^^'B(<^'P) 



xFo,(Ki,^«.-Ai,---,-Ar,a+P,-^^^,---,-^^;^,^,^J. 

fr) 

The second equality follows by the definition of extended Lauricella hypergeometric function j^^^j. □ 



Theorem 3.10. 

pOO pOO T r 

" ^Pj-l„-ti_T:(<i) ;xitl +---+Xrtr;b,d dtl • • • dtr 



^ ;Xltl H hXrtr'.b, d 



j=1 

= r(|3i)---r(|3,)F[;\^^)(a,Pi,--- ,p,;r,xi,--- ,x,;b,d). (75) 
(fH(Pi),--- ,rH((3,) >0] 

Proof 

The demonstration of this theorem depends on the application of the following series identity 
„ mi!---mr! m! 

mi,---,mr=0 m=u 
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specifically, we have 

/ ■■•/ TTtf'^^e-tiiFi'^' (a;r,xiti ++x,t,;b,d]dt, •••dt, 
Jo Jo 

" ^ B(a,y-a) ^^^^ , . . . ^^,1 1 J io ' ' 



mi ,mT=0 j = l 

f r^R ^ ^ r,a ^ , gj^^' (a + mi + ■ ■ ■ + m,,y - a) x^' ---x^-r 

> +m^)■■■^[^r + ^nr]^^ —. ■ 

^ — B I a,y — a 



^b!d (c^ + TTVi +---+mr,y-a) x^' ■■•x^^ 
B(a,y-a) mi!---mr!' 

mi ,mr— 



□ 

In accordance with the previous method, we can obtain the following generalization for classical Lauri- 
cella's hypergeometric function: (a, 3i , • • • , Priyi , • • ■'YT','X■^ > • • • ^'X-r)- 

[v] 

Definition 3.3. Lauricella hypergeometric function (a, (3i , • • • , (3r;yi , ■ • •yr',''! > • • • ^'X-r) is defined by 

^a\ki) ("^iPi,--- >PT;yi,---yT;xi,--- ,xr;b,d) 

mi,X^o B(P,,yi-P,) ^^^^ 

(|xi| + --- + |x,.| < 1; niin{$n(b),$R(d)}> 0) 
When b = d = 0, function (77) reduces to usual Lauricella hypergeometric function fl6| : 

(a,P„...,P.;y„...y.;x,,..,x.]= ^ i-U.-.^. (c,U,...(c.U m.-.-.m.!- 

mi ,---mr— 

Comparing the coefficients of j^^j with those of F)^ , associated with Definition 2.1, it is clear that 
the construction of the coefficients of function (77) corresponds with what we generalize hypergeometric 
function. 

For this new function, we have the following multidimensional integral representation. 
Theorem 3.11. For Lauricella hypergeometric function F^ (a, (3i, • • • , (3r;yi, ■ • -yriX], • • • ,Xr), we have 

r 

Fa'(kO (a,|3i,--- ,PT;Yi,---yT;xi,--- ,Xr;b,d) =]^B((3j,yj -|3j) 

-1 r1 T 



•••^ ri^f'"'n-Uj)^'-P*-'e(^Ki;-^-y-^)(1-XiUl X,U,)-'^dUT.-du,. 

(78) 

(lH(|3j)>0, $n(yj-Pj)>0, j = ,r; min{lR (b) (d)} > 0) 

Proof 

This formula can be easily established by expanding the factor (1 — xiUi — • • • — XrUr) " as 



[1 -XiUi XrUr] " = ^ 



oo f ^ml / \ 

XiUi ) ' • • • XrUr) 

a 



m-i ,TrLr— 



^mi+-+m, rai!---Tar! 



and then integrating out uj, j = 1, • • • ,r with the help of our extended beta integral (11). □ 
The following theorem show that extended Lauricella hypergeometric function F^ j^^^j can be expressed 
as a single integral whose integrand is a product of several extended Kummer hypergeometric functions 
iF5^^'((3;y;z;b,d). ^3 



Theorem 3.12. 



A,(Kl) 



(a,3i,--- ,(3r;Yi,---yT;xi,--- ,Xr;b,d) 



1 



r(a) Jo 



Pi 



;xit;b, d 



-(Kl) 
1 



(3r 
Tr 



;Xrt;b, d 



dt. 



(79) 



Proof 



Taking the series form of iF^'^'' ((3j;yj; Xjt; b, d) , j = 1, • • • ,t and then integrating out t. 



By applying the relation (a] 



mi H hTRr — 1 +mr 



rni H l-m.T_i 



□ 

(a + mi + h mr-i )„ , we can write 



^a'(ki) («^>3i>--- ,PT;Yi,---yT;xi,--- ,Xr;b,d) 



'(|3j,Yj-(3j) 
(Kl) / oc + mi 



mi ,---mr — 1 =0 



i=i 



B(|3j,Yj-Pj) mi!---mr-i 



X2F 



mr-l,|3r 



Yr 



;x,;b,d . (80) 



Theorem 3.13. The Mellin-Barnes type integral representation o/F^ ^^^^j is 
^A^K,) (a, |3i,--- ,(3r;Yi,---Yr;xi,--- ,Xr;b,d) 



1_ /• _ /• r(a- si sr] ^ (|3j -Sj.Yi - Pi) risj] 



(2m)" A, A. Ha) n B((3j,Yj-Pi) 

X (-X1 y ■ ■ ■ i-Xrr'' dSi • ■ • dSr. (81) 

where Li , • • • , Lr are Barnes paths of integration. 
Proof 

The basic ideas are the same when proving this theorem and Theorem 3.4. We take multidimensional 
Mellin transform of function j (a, (3i , • • • , Pr', Yi > ' • ' Yr; — xi , • • • , — Xr; b, d), get the result on the basis 

of applying the method of Brackets and take the inverse Mellin transform of the function we get. 

To obtain the multidimensional Mellin transform of F^'j^^^ j (oc, |3i , • • ■ , Pr', Yi > • ' • Yri ^ ''i > • • • > —'X-r', b, d), 

we multiply both sides of (77) by x^' • • • x^'^^ and integrate with respect to xi , • • • , Xr over [0, oo) x • • • x 
[0, oo). One has 

oo 

T^V^^ ■■■^r'^^'^AdKi] (a,(3l,--- ,Pr;Yl,---YT;-Xl,--- , -Xr', b, d) dxi • • ■ dXr 

r(a + mi +---+m,) -j^ 4^^'((3i+mi,Yi^Pi)("1)"^' ■•■(^1)"^^ 



m I ,■ ■ -TrLT =0 



r(a) J^l B(|3j,Yj-(3j) mi!---m,! 



OO pOO 



= Y_ f (mi, • • • ,mr) 4)m, ■ ■ • Ct'm, (Sl +mi) • • ■ (Sr + mr) . 

mi ■ -mr —0 

The result here is just more comphcated in form than that of equation (56). And with the basic ideas when 
we prove Theorem 3.4 and with the applying of Rule III[1], we can get 

oo /"OO 

X]'"^ ••■x^'-^^f'T' , (a,Pi,--- ,(3r;Yi,---Yr;-xi,--- ,-Xr;b,d)dxi • • • dx^ 







r(a-si s,) ((3j - Sj,Yj - (3j) r(si] 

r(a) y B(pj,Yj-(3j) 
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Taking multidimensional inverse Mellin transform we have 

^a,\ki) («-,|3i,--- ,Pr;Yir--Yr;-xi,--- ,-Xr;b,d) 

^ 1 f f r(cx-si TT 

'[InVrju'Ju Hex) \_\ 

x7*^ • ■ •xr^'^dsi • • ■ ds 



r (CX - Si S,) -py gj^^' (|3j - Sj , Y) - |3j ) r (Sj ) 

L, Ju Hex) \} B(|3j,yj-|3j) 



Theorem 3.13 follows from replacing — xt, i = 1 , • • • , r with xi. 

The way of proving the following two theorems are exactly the same. 

Theorem 3.14. The Mellin-Bames type integral representation ofF^^^^^ is 
Fd!(ki) («:,|3i,--- ,Pr;r,xi,--- ,Xr;b,d) = ^ 



(2m)- 

• /Tr(Pj-Sj)r(sj)B|,^^'((x-si s„T-«) 



r(|3j) B(a,Y-a) 



(-Xl)"'' •••(-Xr)"'^ dSl •••dSr 



w/iere Li , • • • , Lr are Barnes paths of integration. 
Theorem 3.15. For p = q + 1, we have 
q 



r(cxon^^^^^^(fp^q+in"'(p; ::: ;-(xi+-+x.);b,d) 

/ / '^^'^^"^^ Sr)]^B{,^^^ (aj+l -Si Sr, Pj-aj+l) 



1 



(2m 

X r (si ] • • • r (sr) Xi"'^ • • • x^^'- dsi • • • ds 

For p = q 



A r(aj)r(Pj-aj) / ai ... cxq ._ ' 

n r(Pj) V Pi Pq ' +''rj,b,d 

[ ■■■ [ iKd s„P)-a)) 

'•'•J ^Lr ^Lt j^l 



(2m. 

X r(si)...r(sr)x7*' ...x^^'-dsi ..-ds 



5 = 
-Si 



For p < q 



r(Pi)---r(p,)y r(p.+o ' ^ Ui •■• Pq ' ^ ' +^r],b,d 

X r (Sl ) . ■ . r (Sr) X7*' • . • X^*" dSl • . ■ dSr 

where Li , . . . , Lr are Barnes paths of integration. 
Remzirk 

Note that when parameters b = d = 0, integral (83), (84) and (85) reduce to 



nr=ir(aj) 
n?=ir(p,)^''' 



ti)" A, A. n?=i 



■■■ ;-(xi+.-. + x.) 

pi ••■ Pq 

njLl r(«j +Sl + ... + sr 



(2m)'^A, A. nr=ir(Pj+si+... + s, 
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■r(-Sl)...r(-Sr)Xi' ...X^'-dSl ...dSr 



where the contours are Barnes type with indentations, if necessary, such that the poles of F (a, + si + -- - + Sr), j = 
1 , • • • , p are separated from those of V (— sj ) , j = 1 , • • • , r. And this muhiple integral representation for the 
generalized hypergeometric series is given by Saigo and Saxena ^16. . p. 212). <[) 



4. Applications in Hilbert-Hardy type Inequalities 



We have introduced many different generalizations of usual hypergeometric functions, such as, extended 
generalized hypergeometric function, Appell's and Lauricella's hypergeometric functions. A lot of important 
results are established in the previous sections. Now we hope to find some applications in other branches 
of mathematics. In this section, we will use extended Gauss hypergeometric function to establish several 
important inequalities. 

Let f (x) ,g (y) > 0, f (x) e LP (0,+oo) ,g (x) e L^ (0,+oo), 1 + 1 = 1 and p > 1. Then we have the 
following two equivalent inequalities as 

fWg(y) 



OO POO 



y 



dxdy < 



7T 



sm I - 



fP fx] dx 



(87) 



°°lWdxyd,< 



n 



sm I - 



fP fx] dx 



(88) 



where the constant factor 



and 



1 P 



are the best possible in (87) and (88) respectively. And the 



equality in (87) and (88) holds iff f (x) = or g (x) = 0. Inequality (87) is called Hardy-Hilbert integral 
inequality (see [l3|)- During the last decade inequality (87) and (88) were generalized in many different 
ways, some of the generalizations rely heavily on the integral identities involving hypergeometric functions 
and beta functions (see [3|)- [l3|> author obtain some new Hardy-Hilbert type inequalities with 
fractional kernels, which involve the constants expressed in terms of Gauss hypergeometric functions. So it 
is natural to consider whether these new inequalities can be extended to the more general case by using our 
extended Gauss hypergeometric functions. The answer is positive. In the whole process of proving, we can 
get some previously unimagined properties. Note that in what follows, we chose function: O (ki,; z) = exp (z). 
This expression are easy to handle and the final form of our result will thus be more simple. 
We need the following two integral identities. 

Lemma 2. Suppose a, b, c, oc, y, p, q G M are s.t. p^^O, q>0, Q-|-c>b>0 and < a < 2y. Then 



(1 +yxr[]+ccxy 
x'^-i 

(1 +ax)"(1 +yx)' 



exp — yqx 



■ exp 



-aqx 



py 1 



pa 



dx = eP+ V^^B (b, c + a - b) 2F1 
dx = eP+ V^^B (b, c + a - b) 2F1 



a, b y — a _ , 



c,b a — y _ _ 



(89) 
(90) 



Proof 

We start with the integral representation (13) with (ki;z) = exp (z), i.e. 



(1 -t) 



exp I 



dt = B(b,c-b)2Fi 



a, b _ _ 

^ ;z;p,q. 



(91) 



By using the substitutions 1 — t = in (91), we have 



u 



b-1 



1 +u]''"''(1 + u-zu)"''exp -uq--t- du = eP+^B (b, c - b] 2F 



Q, b 



;z;p, q 



(92) 



Then we substitute yx for u in Eq (92) 



(1 +yx)'^ " (1 + (y — zy)x) 



- exp -yqx 



py ^ 

X 
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^ du = eP+^y-''B(b,c-b)2Fi 



Q, b 



;z;p,q 



(93) 



The first assertion of lemma 2 follows from setting y (1 — z) = a and replace c — a with c in above equation. 
To prove the second assertion we can make transformation u = ax in Eq. (92) and then set a (1 — z) —y. 
Replacing c — a with c in the result and then interchange a and c we obtain the relation (90). □ 
Remark 

If we let q = p = 0, lemma 2 reduce to [14, LEMMA 1.] 



dx = y-^B(b,a + c-b)F 



(1 +ax)^1 +Yx)' 



a,b y — a 



a + c Y 



In fact, introducing the factor exp (— ^- — "f^) weakens the symmetry of integrand. So we need two integral 
identities to establish our results. <0> 
It is time to finish the proof of Theorem 2.5. 

The Proof of Theorem 2.5: With the help of the formulas (89) and (90) of Lemma 2 we are able 
to finish the proof of the Euler transformation mentioned in Theorem 2.5. Let us set (\ — and p — ^ in 
the second assertion (90). Then (90) becomes 



exp — yqx 



py- 



dx = e'^+^a^'^B (b,c + a-b)2Fi 



(1 +ax)"(1 +yx) 
By equalizing the identities (94) and (89) we obtain the relation 



c,b _ a — y ap yq 
c + a ' cc ' y ' a 
(94) - 



p-(i-^)q fy 



2r^ 



c,b .1 T.^ Tq 

c + a ay a 



c + a 



;p.q 



(95) 



If we write ^ = 1 — z, then (93) becomes 



c — a,b p ,^ , _ 



c ' '1-z 
By using the notation occured in section 2. we have: 



2r^ 



a, b 



c '1-z 



;p,q 



(96) 



(3i - ai,a2 b 

\z\z ,(1 -z)d 



(3i 



1 -z' 



2ri 



ai,a2 



Pi '1-z 



;b,d 



Substituing equation (95) into (24) we get: 



(97) 



2ri 



ai,a2 



P 



;z;b,d 



:i -zi 



2ri 



ai , P 1 - a2 . -z 



CX2 



' 1 -z 



;d,b 



(l-z)b-zd^l _^iPi-a2-ai 



Pi-a, P,-a2 b 
Pi I — z 



(98) 



This is in fact a generalized Euler transformation. 

Let us come back to the proof of inequalities. Before we state and prove our new inequality, we need to 
define some weighted functions. We define F : (0, oo) ^ R by: 



-q'Ai 



[x + (x^y} lx+ a2y 



„2 exp -aiq dy 

' ^ ' X ai y ' 



where p > 0, q > 0, |j. > 0, ai , a2 > 0, j < < 2, si + S2 > and Az £ 
G : (0,oo) ^ M by: 



1 — S I — S2 _1_ 

q' ' q 



(99) 



^j]- Define 



G (ij;p,q) 



r I ^SI f , iS2 Sxp 

(^ + oi^y) (x + a2yj 



oczy a2 X 



dx 



1 — S I — S2 1 
P' ' P' 



(100) 



By virtue of 



where p > 0, q > 0, fo. > 0, ai , a2 > 0, j < -f^ < 2, si + S2 > and A] e 
lemma 1 we can compute F (x) and G [y) and express them as extended Gauss hypergeometric functions. 
Lemma 3. Suppose |j, > 0, ai , a2 > 0, -j < < 2, ^j, + S2 > 0, Further, let Ai and Az be real parameters 



such that Ai e , and Ai G ^7-^) • If the functions F and G are defined by (99) 

(100) respectively, then 

F(x;p,q)=e'T^k(F; p,q)x^^T^-^^ (101) 



G (y;p,q) = eV^k(G; p,q)y ^ 



(102) 



where 



k(F; p,q) = ai (1 - q'Ai.s, + S2 + q'Ai - 1] <^ 2F 



S2, 1 — q'Ai ai - ai _ _ 
si + S2 ai 



, (103) 



k(G; p,q) = «! " (k^^ (1 -p'Ai,Si + S2 + p'A, - 1 ] <^ 2F1 



Proof 

This lemma is a direct consequence of lemma 2. 
Now we are prepared to derive our main result. 



si, 1 -p'Ai ai - a2 . 
Si + S2 ai 



(104) 

□ 



Theorem 4.1. Let p and q 6e the real parameters such that p>1, q>l, :^ + -^> 1,anrf let p' and q' 
respectively he their conjugate exponents, that is, ^ + ^ = '\ and + = 1 • Define A = :p + and note 
that < A < 1 , for all p and q. Further, suppose > 0, ai , a2 > 0, ^ < -^j < 2, Si + S2 > 0. iff anrf g are 
non-negative measurable functions on (0, 00], then the following inequalities: 



00 /.OO 



exp 



(ai q + a2p) 



+ 



1 X 



^0 



(x + aiy)^'^ (x + a2y) 



^S2 



f (x]g(y)dxdy 



^J2^(l-Sl-S2)+p(Al-A2) 



fP fx] dx 



Jo 



A-(l-Si-S2) + q(A2-Ai ) q 



(y) dy 



and 



(105) 



/ 

Jo 



(si+S2-l) + q'(Ai-A2) 



r 



exp 



-(aiq + a2p) 



1 X 



X ai a2 -y 



(x + aiy]'^' (x + a2y)^^^ 



f (x) dx dy 



XI 



iV(l-Si-S2)+p(Ai-A2) 



fT (x) dx 



(106) 



/loZd for any Ai e (^ ^~'^r'% ^) a?irf A2 G ( ^~^qr^% ^) • T/ie constant 

Kp.q =k(F; q'p,q'q)k(G; p'p,p'q). 

Inequalities (105) and (106) are equivalent. In addition, the equality in (105) and (106) holds iff at least 
one of the functions i or g is equal to zero. 

Proof 

The left-hand side of (105) can be written as 



CO rco 



exp 



-(aiq + a2p)(^ + ^t) 



'0 ^0 
where 



(x + aiy]^" (x+ a2y]^'^ 



f (x)g(y)dxdy 



00 rco 1 



^0 



Ki"' (x,y)Kf (x,y)K^-^(x,y)dxdy 

(107) 



FT-q' (x;q'p,q'q)fP (x) xi"^^ 



(x+aiy)'^ (x + a2y)'" y"'^ 



'x 



yx- 
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K2 (x,y) 



I , ISW , .S2 rj'A. '^^P ~'^2 P q r _1 



X , pai 
--P — 



K3 (x,y) = xP'^' FP (x; q'p, q'q)y^^^G^ (y;p'p,p'q) W g'' [y] . 

(It is easy but time-consuming to verify Eq (107).) and the functions F and G are defined by (99), (100). 
Now, since the exponents satisfy identity: ■^ + -^ + {1— \) — ^■ We can use Holder's inequality, namely, 

00 nOO ] ] 

/ {x,y)Kf (x,y)Ki-^(x,y)dxdy 

Jo 

<l I / K, (x,y)dxdy / / K2(x,y)dxdy / / K3(x,ij)dxdy 

Jo J \Jo Jo / \Jo Jo 

1 



< 



xP-^iFP (x;q'p,q'q)fP (x) dx 



y^-^^GMyip'p.p'qjgMyldy 



l-A 



< 



^"yi'^^G'' (ij;p'p,p'q) [y] dyj ^ ^^^^'^^ (x; q'p, q'q) fP (x) dx 
|eT^{q'P+q'^)kP (F;q'p,q'q)}^ {el^^f^'P+P'^^k^ (G;p'p,p'q)}^ 



^JV(l-s,-s.)+p(A,-A,)^p, 



where 

{ef^('''P+''''^)kP{F;q'p,q'q)}^{e^(p'P+P''^^kMG;p'p,p'q)}^=e2(P+'i'k(F;q'^^ 

It marks the end of the proving of the first part of the theorem. The equivalence of (105) and (106) and 
the conditions for equahty to hold can be proved in the way which explained in |14l p. 649, Theorem 1]. □ 
Remark 

In fact, the way we prove the inequalities (105) and (106) provides a common method to generalize 
inequalities related to gauss hypergeometric functions. Once some integral identities are established, we can 
prove it according to the primary way of proving. For instance, the inequalities in the reference [4] can be 
generalize in the same way. (we should notice that we use definition 2.2 instead of 2.1). 

On another aspect, we can directly use the definition of extended hypergeometric functions to construct 
inequalities. because the new definition is strong and flexible enough to get the desired results. <) 



5. Conclusion 

By using extended Beta functions we have defined some extensions of generalized hypergeometric func- 
tion and Lauricella's hypergeometric function. In some special cases, those extensions have been studied 
in Q, [111 and With flexible manipulating of extended definitions, we have obtained some generaliza- 
tion of properties of hypergeometric functions in classical theories. We have also explained an example of 
Hardy-Hilbert type inequality involving extended Gauss hypergeometric functions in detail, which may be a 
precedent for applying these extension in other branches of mathematics. It is worthwhile to point out that 
some fundamental relations between fractional calculus and extended generalized hypergeometric functions 
have been proposed in this paper. If we want to study it in the aspect of fractional calculus, we need to get 
more properties of the fractional operator defined by (43). As far as this is concerned, we have gotten some 
ideas, which will be presented in next publications. 
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